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The electron-hole symmetry in the structure "graphene - insulating substrate -semiconductor 
gate" is violated due to an asymmetrical drop of potential in the semiconductor gate under positive 
or negative biases. The gate voltage dependencies of concentration and conductivity are calculated 
for the case of Si02 substrate placed over low- (moderate-) doped p-Si. Similar dependencies of the 
optical conductivity are analyzed for the case of high-K substrates (A1N, AI2O3, HfC>2, and Z1-O2). 
The comparison of our results with experimental data shows a good agreement for both cases. 
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The energy spectrum and scattering mechanisms in 
graphene are symmetric with respect to electron-hole re- 
placement, so that the transport phenomena are identi- 
cal for the electron and hole doping cases, see reviews 1 
and 2. According to the theoretical analysis, as well as 
the transport and optical measurements (see review [3[ 
and references therein), the Coulomb renormalization in 
graphene leads to weak (~10 % for typical parameters) 
asymmetry of electron and hole responses. These results 
were obtained for the electrostatically doped graphene 
structures with metallic (or heavily-doped semiconduc- 
tor) gates. The sheet concentration of carriers, n g , is de- 
termined by the gate voltage V g according to the plane 
capacitor formula 



a(V g - V ), 



47r|e|d' 



(1) 



where e s is the dielectric permittivity of substrate of 
thickness d, see Fig. la, and voltage Vq corresponds to 
the electroneutrallity point; below we suppose Vq — > 0. 
Thus, transport phenomena should be symmetric with 
respect to the sigh flip of V g (electron- hole symmetry). 
But this symmetry appears to be violated for the struc- 
tures with low- (moderate-) doped semiconductor gates. 
It is due to differences in the drop of potential and charge 
distributions for positive or negative biases, when the de- 
pleted or heavily-doped regions appear, see Figs, lb and 
lc, respectively. This simple mechanism was not ana- 
lyzed yet and it is timely to re-examine Eq. (1) for the 
case of low-doped semiconductor gates, which were used 
in several experiments. 041] 

In this letter we solve the electrostatic problem for 
graphene structures with thick Si02 or thin high-K (A1N, 
AI2O3, Hf02, and Zr0 2 ) substrates placed over p-Si of 
different doping levels. The gate voltage dependencies of 
the concentration and the dc conductivity as well as the 
interband optical conductivity are calculated and com- 
parisons with the experimental data [4-6] are performed. 
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FIG. 1: (a) Structure schematic for back-gated graphene 
(grey) with dielectric substrate of thickness d and permittiv- 
ity e s over p-doped Si of permittivity e. (b) and (c) Potential 
(upper) and charge (lower) distributions along semiconductor 
gate for the electron and hole doped Si cases, respectively. 



The distribution of the electrostatic potential cp z across 
the structure under consideration (along z-axis, z < d) 
is governed by the Poisson equation 



difi z 
dz V dz 



< z < d 
z < 



(2) 



where e z stands for the dielectric permittivity and p z is 
the charge density in graphene (at z — > d) and in Si gate 
(at z <0). Using the Gauss theorem at z — > d — and 
connecting the derivatives of ip x at z —> one obtains the 
boundary conditions 



dz 



4tt 



z=d-0 



d(fi z 
dz 



+0 



(3) 



Here <j a 



rd+0 , 

Jd-n dz P* 



=F |e|n s is the charge density in- 



duced in graphene which is written through n g and — or 
+ correspond to the electron or hole doping. The linear 
potential distribution ip z = tp z =o + A-n(j g z/ e s takes place 
in the substrate region d > z > 0, where p z — 0. In 
the gate region z < Eq. (2) is transformed into the 
second-order differential equation 
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which is written for p-Si (n-Si) through the differences 
between acceptor (donor) and hole (electron) concentra- 
tions, ha (fifl) and n^z (n ez ). The boundary conditions 
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FIG. 2: Concentration of carriers n 9 (electrons if V g > or 
holes if V g < 0) versus gate voltage V g for graphene on SiC>2 
substrate of 300 nm or 100 nm width (solid or dotted curves) 



placed over p-Si gate with doping levels: ua 
(1,4), 10 16 cm" 3 (2,5), and 10 17 cm" 3 (3). 
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for Eq. (4) are <p z ^~oo = and the right Eq. (3), which 
determines an electric field at z — > — 0. 

Below, we restrict ourselves by the case of p-doped 
Si, where the non-degenerate hole approximation n^ z = 
ha exp(eip z /T) is valid for ua < 5 x 10 18 cm~ 3 at room 
temperature T =300 K. For higher ua or for the n-doping 
case, nhz or n ez depend weakly on tp z and deviations from 
Eq. (1) are weak. Since Eq. (4) does not depend on z 
explicitly, we obtain the first-order nonlinear equation 



e dip z 
T dz 



2 



ey z 
T 



1 



(5) 



where the characteristic length zt is determined by 
the relation (47re 2 / eT)n a z\ — 1. Further integration 
gives the potential and charge distributions, which are 
schematically shown in Figs, lb and lc, and <p z =o is 
connected with n„ through the boundary condition at 



-0: 



2 \ n T 



e<p z = 
T ' 



(6) 



Here we introduce the characteristic 2D concentration 
n T = n A z T oc \JtiaT which is about 0.8 x 10 10 cm~ 2 
if n A = 10 15 cm" 2 and T =300 K. Neglecting the low- 
doped region n g < nx, we obtain the simple relation 
between V g and n g : 



aV a 



TLj 1 - 



c s zt 



(n g /n T ) 2 /2, (h) 
-21n(n fl />/2nr), (e) ' V ' 



which generalize the standard plane capacitor formula 
given by Eq. (1). For the electron-doping case (V g > 0), 
both oc n g and oc n 2 contributions are essential depend- 
ing on the dimensionless factor e s ZT/ed. By contrast, for 
the hole-doping case (V g < 0) the ln-correction can be 
only detected in the low-doped graphene. Note, that 
there is no temperature dependency for V g < and 
these dependencies are weak (appears in ln-correction) 
for V g > 0. 

Fig. 2 shows the carriers concentrations versus V g in 
graphene placed over the SiC^/p-Si structure depend- 
ing on doping level and substrate thickness, ua and d. 




FIG. 3: (a) Fitting of experimental data from Refs. 4 and 
5 marked as 1 and 2: (a) momentum relaxation rate v v ver- 
sus energy vp and (b) conductivity a versus V g . Squares and 
triangles are experimental points from Refs. 4 and 5, respec- 
tively. 



These dependencies are substantially different for e- and 
/i-doping regimes. If V g < 0, the plane capacitor ap- 



proximation (1), with the n^-independent concentration 
n g oc d , is valid. If Vg > 0, the dependency n g (V g ) 
varies from the square-root to linear function and the 
electron-hole symmetry is restored if ua exceeds 10 17 
cm -3 at d =300 nm. For thiner substrates (d =100 
nm is shown), the symmetry is restored at higher con- 
centrations. The square-root dependency for low-doped 
p-Si appears due to the depleted region of thickness ddopi 
(see Fig. lb), which is determined by the charge neu- 
trality condition, n s ~ n^rfdcpi- Using Eq. (1) with 

>V a for the case 



d — > d + (idepi, one obtains 
ddcpl ^> d. 

The asymmetrical electrostatic doping of graphene un- 
der sign-flip of V g leads to similar changes of the gate 
voltage dependencies of static conductivity, cr(V g ). These 
dependencies in structuresplaced on p-Si with ua — 10 15 
cm~ 3 were measured 0, [f| under investigations of the 
Drude and interband absorption in large-area samples. 
In order to fit cr(V g ) we use the phenomenological mo- 
mentum relaxation rates v p suggested in Ref. 7 with the 
parameters taken from the experimental data of Refs. 4 
and 5 for the hole doping cases V g < 0, when Eq. (1) 
is valid, see Fig. 2. Using the relaxation rates shown in 
Fig. 3a and the dependencies n{V g ) for V g > 0, we obtain 
an excellent agreement between u{V g ) and the results of 
Refs. 4 and 5, as it is shown in Fig. 3b. The asymmetry 
of a(V g ) is more pronounced than n(V g ) (c.f. Figs. 2 and 
3b) because of an additional energy dependency of v p . 

Further, we turn to consideration of the structure with 
a thin high-K substrates placed over p-Si. For these struc- 
tures, Eq. (7) contains the parameter £ s zt jtd 3> 1 and a 
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FIG. 4: (a) Concentration n g versus V g for graphene on A1N 
(black), AI2O3 (red), Hf0 2 (green), and ZrO a (blue) sub- 
strates of width 10 nm placed over p-Si of doping level 10 5 
cm -2 , (b) and (c) The same for p-Si of doping levels 10 16 
cm -2 and 10 17 cm -2 , respectively. 



deviation from the standard formula (1) increases. Sim- 
ilar to Fig. 2, the dependencies of concentrations n s on 
V g are shown in Fig. 4 for A1N, A1 2 3 , Hf0 2 , and Zr0 2 
substrates [8j of width 10 nm at different doping levels. 
Here the square-root behavior of n(V g ), which is weakly 
dependent on e s , takes place at V g > 0. For the hole dop- 
ing regime, concentrations n(V g ) deviate from the linear 
functions at low \V g \, due to the ln-contributions in Eq. 
(7), and curves are different for different substrates due to 
variations of e s . For higher doping levels riA, the induced 
concentration n g increases with V g faster, c.f. Figs. 4a- 
4c where ua varies from 10 15 to 10 17 cm -3 . The interval 
of V g applied is restricted by the breakdown condition 
V g /d <5 - 10 MV/cm, [9| so that we consider the region 
\V g \ <5V. 

The essential asymmetry of electron-hole doping 
changes conditions for interband absorption because of 
a different Pauli blocking effect for V g > and < 0. At 
low temperatures (or for heavily-doped graphene) the ab- 
sorption is blocked in the spectral region hu < 2ef- Here 
the Fermi energies of electrons and holes £j? a v/^g °- e " 
pend on V g asymmetrically, see Figs. 2 and 4. So that, 
the threshold of absorption spectra for V g > and < 
should be different: at fixed Kto, higher threshold bias is 
necessary for the electron doping regime, i.e. the asym- 
metric Pauli blocking takes place. In Fig. 5a we plot the 
electron and hole Fermi energies versus the gate voltage 
for structures with AI2O3 substrate of different widths 
placed over p-Si of different doping levels. Similarly to 
Fig. 4, there is a visible dependency on ua if V g > 
while for Vg < the dependency of Ef on d is only essen- 
tial. The dependencies SF(V g ) for structures with the all 




FIG. 5: Electron and hole Fermi energies versus gate voltage 
V g for AI2O3 substrate of widths d =7 and 15 nm (solid and 
dashed curves) at p-Si doping levels 10 15 , 10 16 , 10 17 , 10 18 
cm -3 marked 1 to 4, respectively, (b) The same for different 
substrates [A1N (black), AI2O3 (red), Hf0 2 (green), and Zr0 2 
(blue)] of widths d =10 nm at p-Si doping levels 10 15 and 10 17 
cm -3 (solid and dashed curves, respectively). 



high-K substrates under consideration are shown in Fig. 
5b for low- and heavily doped p-Si. Once again, Ef(V 3 ) 
is not dependent on e s at V g > and the dependencies 
on tia are negligible at V g < 0. Note, that the threshold 
of absorption appears at low voltages for the hole doping 
case, V g > — 1 V if fkv <0.6 eV. For the electron doping 
case, there is no threshold up to the breakdown voltage, 
if huj >0.3 eV. 

An essential distinction between the spectral depen- 
dencies of interband absorption for the electron and hole 
doping regimes was reported in Refs. 4 and 5 for the 
structures of large-area graphene placed on Si0 2 sub- 
strate over low-doped p-i. These results are in a quali- 
tative agreement with the asymmetrical dependencies of 



£F(Vg) which can be obtained form Fig. 2. But a quanti- 
tative description of the dynamic conductivity a u should 
take into account both the temperature broadening and 
the disorder effects, similar the consideration of the hole 
doping regimes performed recently, [l0| and it requires a 
special consideration. Here we consider the asymmetrical 
threshold of interband absorption in the graphene-based 
modulator with the Al 2 03 substrate of width 7 nm and 
the waveguide designed for the telecommunication spec- 
tral range, Tllo ~0.8 eV. [6| The depth of modulation, 
AT U , should be proportional to the carrier- induced con- 
tribution of the dynamic conductivity, AT W cx ReAcr^, 
where ReAa^ is given by [loj 



ReAcr,. 



Cintc 



exp[(hu-2e F )/2T} + l 



G intc 



e 



(8) 



Here we consider the case ef ^> T, so that the chemical 
potential is replaced by Ef- Using £f{V 9 ) shown in Fig. 
5, below we consider the dependency of AT U cx ReAcr w 
on V g normalized to the unit amplitude of modulation. 

The normalized conductivity 2Re<7 w /t7i ntcr is plotted 
in Fig. 6 for the conditions of Ref. 6 at different doping 
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FIG. 6: Dissipative interband conductivity ReAcT^, normal- 
ized to o-jntcr, versus V g for graphene on AI2O3 substrate at 
p-Si doping levels 10 15 (1), 2 x 10 18 (2), and 3 x 10 18 cm" 3 
(3). Triangles show experimental data for AT W from Ref. 6 
which are normalized to the amplitude of modulation. 

levels of p-Si. One can see, that the experimental data for 
the normalized dependency of AT^ on V g can be fitted 
at n A 3 x 10 18 cm -3 and T =300 K (a doping level is 
not given in [6J and we use ua as a single fitting param- 
eter). The thresholds of absorption are in good agree- 
ment with the measurements both for the electron and 
hole doping regimes but the temperature-induced smear- 
ing leads to more abrupt jumps in comparison to the 
experimental spectra. It is because a disorder contribu- 
tion did not taken into account, see similar analysis of 
the interband absorption spectra [1, @ in Ref. 10. In ad- 
dition, the narrow waveguide structure forms a strip-like 
capacitor with the non-uniform charge distribution across 
the waveguide; [ll[ this leads to an additional smearing 
effect. There is a weak doping dependency of AT W at 



V g < and at Vg > the modulation condition Tioj ~ 2ef 
satisfies if ha > 10 18 cm" 3 . The last estimate is valid 
for the in-plane homogeneous structure, while for a bus 
waveguide 6] doping level should be lower. A complete 
description of the results of Ref. 6 requires to take into 
account both the disorder contributions and the waveg- 
uide geometry 

Now we list the assumptions used in our calculations. 
The consideration of the non-degenerate holes is valid for 
any ua and any V g because of formation of the depiction 
region at V g > and weak corrections to Eq. (1) at V g < 
0. Apart from the mechanism under consideration, the 
Coulomb renormalization in n- or p-doped graphene leads 
to an additional (up to 10 %) asymmetry of response. [3, 
[l2} Such a contribution can be neglected in the structures 
analyzed here (with a low-doped p-Si gate or a thin high- 
k substrate). But in typical samples (with the electron- 
hole asymmetry <10%) both the Coulomb effect and the 
electrostatic contribution should be taken into account. 
This could result in decrease or increase of the electron- 
hole asymmetry. In the vicinity of the zero bias, V g ~ 
0, the long-range disorder [l3[ should modify the effect 
under consideration. 

Summarizing, the mechanism of the electron-hole sym- 
metry breakdown, caused by different drops of potentials 
in low-doped Si for opposite signs of bias, have been an- 
alyzed and the results are in agreement with the recent 
transport and optical measurements. These results are 
also important for device applications such as the field 
effect transistors [l4[ or the high-frequency multipliers. 
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